Phase separation in a chaotic flow 
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The phase separation between two immiscible liquids advected by a bidimensional velocity field is 
investigated numerically by solving the corresponding Cahn-Hilliard equation. We study how the 
spinodal decomposition process depends on the presence — or absence — of Lagrangian chaos. A 
fully chaotic flow, in particular, limits the growth of domains and for unequal volume fractions of the 
liquids, a characteristic exponential distribution of droplet sizes is obtained. The limiting domain 
size results from a balance between chaotic mixing and spinodal decomposition, measured in terms 
of Lyapunov exponent and diffusivity constant, respectively. 
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A system of two immiscible fluids at rest will gradu- 
ally phase-separate, forming domains whose size grows 
algebraically with time. Everyday experience, however, 
shows that by continuously stirring or shaking the fluids 
the domains or droplets of the phases (say, oil and vine- 
gar) break and coalesce, leading to a dynamic stationary 
state with domains of finite size. 

A first approach consists in modelling this situation by 
subjecting the binary fluid to a homogeneous shear veloc- 
ity field [IJ. However, experiments ||], numerical simu- 
lations |3], and more recently analytical approaches 
show that in such a situation infinitely long domains 
aligned with the flow are formed. The effect of the ve- 
locity field is to counter the Rayleigh instability, stabil- 
ising lamellar and (in certain cases) even cylindrical do- 
mains [||J^]. Domain breakup in those situations takes 
place only at large Reynolds numbers, and is generally 
attributed to inertial effects [jl|,^| . Studying these inertial 
effects numerically is difficult, as a realistic description of 
the feedback of domain shape on the flow is required |Q . 
The corresponding simulations are therefore limited by 
finite size effects. 

In this paper we investigate a different mechanism by 
which domains of finite size can be stabilized in a demix- 
ing system. In particular, we show that a saturation of 
the average length scale takes place even in the absence 
of inertial effects if the flow has Lagrangian chaos, (i.e. 
if the trajectories of nearby starting points diverge expo- 
nentially with time). This is interesting for two reasons: 
First, with an appropriate time-dependence of the veloc- 
ity field one can still have Lagrangian chaos in a situation 
of high viscosity in which inertial effects are negligible — 
this is how one mixes pastes. Secondly, it is possible in 
that case to decouple the hydrodynamic problem from 
the phase separation. This problem of a passive, phase 
separating scalar field is of course much simpler, so that 
simulations using large systems are possible. Our ap- 
proach therefore extends earlier extensive studies of pas- 
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sive scalar advection by periodically driven chaotic flows 

§■ 

Our study is also related in spirit to earlier studies of 
advection by 'synthetic' velocity fields tuned to model 
turbulent flows [§,0. Phase separation was studied in 
this context in Ref. An essential difference, com- 

pared to our work, is that in such turbulent flows the 
separation between nearby tracer particles appears to in- 
crease algebraically, rather than exponentially, with time. 

We consider a two dimensional flow that can be tuned 
to be regular, mixed, or fully chaotic. Specifically, the in- 
compressible velocity field v(x, y, t) is a modified version 
of the so-called time-dependent Harper map Q (related 
to the 'partitioned-pipe mixer', a special case of 'egg- 
beater flow' [||). The dynamics takes place on a square 
of side L, with periodic boundary conditions. The veloc- 
ity field is an alternating sequence of shears in the x and 
in the y direction with a time period T, 
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The parameter a controls the chaoticity of the trajecto- 
ries. If a is small, the two semicycles are composed into 
the smooth, laminar velocity field: v x = — 1 ^ L sin (^J^); 
v y = ^j^- sin (^p). For larger values of a the trajecto- 
ries stretch and fold, and the flow becomes chaotic. In 
order to visualise this, it is convenient to follow the po- 
sition of a point at the end of each cycle. This 'kicked 
Harper' map is shown in Fig. |l| for several values of a. 
For a ~ 0.2 the flow is a mixture of laminar and chaotic 
regions, and becomes fully chaotic around a ~ 0.4. In 
the chaotic situation, it is convenient to characterise the 
flow by the Lyapunov exponent A, defined by the fact 
that nearby starting points separate as ~ e xt . We have 
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computed A as in Ref. and found that the relation 
A c± 1.96 ln(3.35a)/T is a good approximation through- 
out the chaotic regime, a > 0.4. 
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FIG. 1. The maps obtained from snapshots at intervals T 
of the lines of current with the dynamics ([!]), starting from 
various initial conditions. Figures for a = 0.1 (top left), 0.25, 
0.4 and 1.0 (bottom right). 



The spinodal decomposition of the two-component 
fluid is described by the Cahn-Hilliard equation 



d<t>{r,t) 
dt 
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Here (j) is a dimensionless concentration field, the concen- 
trations of the species are [1 ± 4>]/2. We work at T = 0, 
since temperature is irrelevant in this process ^3[ . The 
free energy functional is of the Ginzburg-Landau form 
and reads 



F[4>] = I 
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Here, £ is the equilibrium correlation length controlling 
the width of the interfaces, and d is the number of spatial 
dimensions. We consider two topologically different sit- 
uations (j) {(j)) 7^ : a species is less abundant than the 
other and forms disconnected droplets, and (it) ((/)) = 0: 
the two phases are in equal quantity and form a bicon- 
tinuous structure. Situation (%) has been studied experi- 
mentally @. 

In a chaotic flow, the passive scalar mixes rapidly, 
whereas in the case of phase separation this tendency 
is opposed by surface tension. The competition between 
these two effects can be quantified through two adimen- 
sional parameters, D = TT/£ 2 (the adimensional trans- 
port coefficient of the Cahn-Hilliard equation) and the 
chaoticity parameter a, or alternatively the adimensional 



Lyapounov exponent AT. A large D means that appre- 
ciable diffusive transport will take place during each lam- 
inar half cycle. A large AT, on the other hand, means 
that the mixing process is efficient within a few cycles. 
Note that A can also be interpreted as an average elon- 
gation or shear rate experienced by the fluid particles. 

Equation (^J) is integrated numerically with the veloc- 
ity field ([j]), using the implicit spectral method developed 
and discussed in Ref. |H| . The results are presented with 
time and length units chosen as the cycle period T and 
the interfacial thickness £, respectively. The system size, 
lattice parameter and time step are L = 512£, Ax = £, 
and At = 5 • 10 _4 T respectively. 






FIG. 2. Evolution of an assembly of large droplets in the 
chaotic flow for times t = 0, 0.65, 1.50 and 11.8. 



Existence of a stationary state. — We first show that a 
purely chaotic flow does indeed stop the domain growth. 
In Fig. [| we show the evolution of a phase-separated 
sample with ((/>) = 1/2, upon turning on a chaotic ve- 
locity field (a = 0.4). The large droplets of the initial 
configuration are broken into smaller droplets, until a sta- 
tionary state where droplets successively grow and break 
is reached. Fig. || shows the late stages of coarsening of a 
system with equal concentrations of phases (((j)) = 0) in 
the four velocity fields of Fig. |l|. For a — 0.1, the velocity 
field is laminar. We observe in that case structures very 
similar to those found in a homogeneous shear flow, but 
which now follow the winding flow lines. In the mixed 
case, a = 0.25, large-domain structures form in the lam- 
inar regions of the flow, and break into very small do- 
mains in the chaotic ones. In the fully chaotic situation, 
a > 0.4, a dynamical stationary state is reached, with 
small domains continuously breaking and reforming. For 
a = 1.0, the sinusoidal nature of the underlying velocity 
field becomes apparent. This snapshot nicely illustrates 
the typical 'stretch and fold' processes characteristic of 
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chaotic advection [||. 




FIG. 3. Late stages of coarsening in the velocity fields of 
Figure [lj in the same order. 

Scaling properties in the chaotic flow. — In the sta- 
tionary regime, it is clear from Figs. || and || that there ex- 
ist a typical length scale L* which depends on the param- 
eters D and A: this will be confirmed below by a quantita- 
tive analysis. As in the pure coarsening case |l3|] , scaling 
properties are expected in the regime £ <C L*(D, A) <C L. 
The length scale L* may be estimated by the following 
simple argument. In the absence of flow, the domains 
grow as L(t)/£ ~ (Dt/T) 1 / 3 , and this growth is stopped 
by the chaotic flow which introduces a time scale A -1 . 
Hence, we estimate L* ~ L(t = A -1 ) and predict 

/ D\ 1/3 

L*(D,\)~t[—) . (4) 

Isolated droplets, (0) ^ 0. — Following Ref. [l4|| , we 
characterize the assembly of droplets by computing the 
distribution of droplet surfaces f(S), where f(S)dS is 
the probability that the surface occupied by a droplet is 
between S and S + dS. In the scaling regime, we expect 
this distribution to be of the form 

where S* is a typical droplet area. In Fig. |], data 
obtained for a wide range for the values of D and A 
are collapsed by using a reduced variable S/S*, with 
S* (D/X) 0,62 . The data collapse is satisfactory, and 
the result for S* reasonably close to what would be ex- 
pected from Eq. (f|), i.e. S* ~ (D/\) 2 / 3 . Finding an 
exponent slightly smaller than the one expected theoret- 
ically is not surprising, since the typical domain sizes are 
rather small (S* < 50£ 2 ), so that the asymptotic value 
for the domain growth exponent in the absence of flow 



may not be reached. The rescaled distribution functions 
exhibit an exponential tail, J~(y) ~ e~ y (dotted-dashed 
line in Fig. Such distributions are very is similar 
to those found in the experiments of Ref. [Q. For the 
largest droplets, deviations from the exponential fit are 
observed, indicating either insufficient statistics or a dif- 
ferent scaling behaviour for the extreme values of S . 
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FIG. 4. Surface distribution of the droplets rescaled ac- 
cording to Eq. (||), where the typical surface S* is given by 
S* ~ (D/X) ' 62 . The dotted-dashed line is a fit to an expo- 
nential form, and the data are for a range a G [0.4, 3.0] and 
D G [100, 2000]. 
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FIG. 5. Log-Log plot of the typical domain size L* as a 
function of the ratio D/X for the same values for D and a as 
in Fig. |. The full line has a slope 0.27. 

Equal concentrations: (</>) = 0. — In the case of 
equal concentrations, the domains are ramified and ex- 
tend throughout the sample, so that the area is not a use- 
ful measure of domain size. A characteristic domain size 
can nevertheless be obtained from the two-point corre- 
lation function C(r,t) = L~ 2 J d 2 x((j)(x, t)<j)(x + r, t)), 
which is the Fourier transform of the structure factor 
measured in light scattering experiments. Performing a 
time average over many configurations shows that the 
bicontinuous structure is on average perfectly isotropic, 
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as it is in the absence of flow. One can therefore av- 
erage C(r,t) over orientations to obtain a one variable 
function, C(r). The characteristic domain size L* can be 
defined by C(L*) = 0.5. Fig. || displays this domain size 
for various combinations of D and A, as a function of the 
ratio D/X. The data can be fitted by L* ~ f (D/XT) ' 27 
Again, this is in reasonable agreement with the scaling 
analysis, Eq. (|]). Larger simulations, with smaller val- 
ues of A, would be necessary to obtain larger domain sizes 
and avoid the crossover effects which are well known in 
spinodal decomposition simulations p~3|| . 




3 

r/D - 27 

FIG. 6. Two-point correlation function C(r) at fixed a for 
for various mobility D, as a function of the rescaled variable 
r/D ' 27 . MAIN: a = 1.0 and D = 100, 200, 400, 1000 and 
2000. INSET: a = 0.4 and D = 100, 200 and 400. 



A more detailed information on the domain structure 
is obtained from the full correlation function C(r). Here, 
one could expect from the scaling hypothesis a behaviour 
of the form 



C(r)*c(£) 



(6) 



with C a universal function. This hypothesis is tested in 
Fig. |^, where C(r) is represented for fixed A and vari- 
ous values of D. At fixed A, a good collapse of the data 
obtained for different D is achieved by using a rescaled 
variable r/D - 27 . The inset of Fig. ||, however, shows 
that the shape of the scaling function slightly depends 
on A, so that the universal scaling expressed by Eq. (||) is 
not valid. We attribute this change of the scaling func- 
tion with the flow pattern to the fact that even in the 
chaotic regime the flow cannot be considered as being 
homogeneous and isotropic, but exhibits an underlying 
sinusoidal structure. This is in contrast with the case of 
isolated droplets where the droplet distribution was not 
affected by this structure, recall Fig. Q 

We have studied the phase separation in conditions in 
which the species boundaries are passively advected by a 
incompressible flow. We have shown that a chaotic flow 



results in a steady state with domains of finite size re- 
sulting from the balance between spinodal decomposition 
and chaotic advection, Eq. (||). This should be contrasted 
with the situation observed in turbulent flow, where the 
flow intensity must exceed a threshold in order to stop 
domain growth . Such a difference can be traced back 
to the fact that Lyapounov exponents for passive scalar 
advection are actually in the latter case. The essential 
approximation in our work, compared to realistic experi- 
mental situations, is the assumption that the flow pattern 
is not modified by the domain growth. This assumption, 
however, may not be unrealistic if the two fluids have 
similar viscosities and if the capillary stresses are small 
compared to viscous stresses. This is measured by the 
capillary number C a = rjX/fa/L*), where r\ is the viscos- 
ity and 7 the surface tension. In highly viscous fluids, 
C a is expected to be large, so that the decoupling is pos- 
sible. This decoupling also makes it possible to consider 
analytical treatments. 

We acknowledge useful discussions with A. J. Bray, B. 
Cabane, P. Leboeuf, J. F. Pinton and J. E. Wesfreid. 
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